Mechanical squeezing via parametric amplification and weak measurement 
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Nonlinear forces allow motion of a mechanical oscillator to be squeezed below the zero-point 
motion. Of existing methods, mechanical parametric amplification is relatively accessible, but pre- 
viously thought to be limited to 3dB of squeezing in the steady state. We consider the effect of 
applying continuous weak measurement and feedback to this system. If the parametric drive is 
optimally detuned from resonance, correlations between the quadratures of motion allow unlimited 
steady-state squeezing. Compared to back-action evasion, we demonstrate that the measurement 
strength, temperature and efficiency requirements for quantum squeezing are significantly relaxed. 



Recent experiments have demonstrated impressive 
progress in cooling towards the ground state and mea- 
suring the zero-point motion of mechanical oscillators. 
This brings within reach the observation of nonclassical 
phonon states, with applications in quantum informa- 
tion and tests of quantum mechanics [li|. The most suc- 
cessful systems to date involve cryogenically cooled high 
frequency oscillators strongly coupled to optical or mi- 
crowave fields [2, 3]. However, techniques to manipulate 
quantum states and investigate nonclassical behaviour of 
phonons, apart from creating single phonon statesf4], are 
less well developed. 

A squeezed state, in which the variance of one quadra- 
ture of motion is below the zero-point motion, is the most 
accessible of quantum resources in optomechanical sys- 
tems. This can be achieved, for example, by resolved 
sideband cooling using squeezed or modulated input 
light [5|, |6|. Also promising is squeezing via back-action 
evading measurement (BAE)[7|, which is close to being 
realized[8|. These schemes would allow for ultra-sensitive 
force detection 9] and normal mode entanglementtlo!| but 
are constrained by the requirement of strong coupling 
to the optical mode. Additional downsides are the re- 
quirement of ultra- low temperatures and the side-effect of 
parametric instability due to strong radiation pressure [8|. 

Electrostatic forces, on the other hand, are strong 
enough to create nonclassical states in an oscillator by 
driving it into the nonlinear regime |1|. It has been pre- 
dicted that micro- and nano-electromechanical systems 
(MEMS/NEMS) can be engineered in this way to ex- 
cite arbitrary Fock states [Uj and induce macroscopic 
quantum tunnelling [l3|. In a similar fashion, mechani- 
cal squeezing can be achieved via mechanical parametric 
amplification (MPA)|13|. MPA exploits nonlinearities in 
the electrostatic driving field [IJ], the resonator's intrin- 
sic motion[15| or a coupled charge qubitjla]. A periodic 
modulation in the oscillator's spring constant at twice 
its resonance frequency gives rise to an in-phase am- 
plified quadrature and an out-of-phase damped quadra- 
ture. The amplified gain approaches infinity at thresh- 
old, whereas the squeezing due to damping is limited to 
a factor of one half. This is a long-standing problem that 



also limits the intracavity variance of an optical para- 
metric oscillator 17]. While there are non-equilibrium 
schemes that allow stronger squeezing in MFA'ISI, |19|, 
to our knowledge there have been no previous proposals 
to overcome this limit in the steady state. 

Here, we consider a mechanical oscillator with de- 
tuned parametric driving and continuous quantum mea- 
surement. Such a scheme is relevant to MEMS/NEMS 
and optomechanical systems subject to electromechani- 
cal forces, as shown in Fig. [TJi. Our approach is sim- 
ilar to that of '201, where resonant parametric driving 
and resolved sideband cooling are combined with quan- 
tum measurement, allowing the inference of mechanical 
squeezing. Our detuned MFA (or DMFA) approach is 
illustrated in Fig. ^p. The parametric drive amplifies the 
noise of one quadrature of motion, increasing the infor- 
mation obtained by continuous weak measurement. De- 
tuning the drive from resonance causes the quadratures 
to rotate such that the measurement of the amplified 
quadrature contains information about the fluctuations 
of the squeezed quadrature, thereby further reducing its 
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FIG. 1. a) Schematic of DMPA for an optomechanical sys- 
tem with capacitive modulation of the spring constant k. The 
measurement is mixed at oJm -I- A to provide the two position 
quadratures, b) Phase space diagrams (not to scale) depict- 
ing unequal variance of orthogonal quadratures as noise el- 
lipses. Shown is a thermal state in a frame rotating at ojm 
as a parametric drive, measurement and detuning are turned 
on. Darker ellipses indicate conditional variance. 



conditional uncertainty. 

We find that while measurement back-action degrades 
the squeezing from a resonant parametric pump, a con- 
tinuous weak measurement enhances the conditional 
squeezing if the pump is optimally detuned. This can 
be made unconditional by applying appropriate negative 
feedback based on the measurement [2 ll| . Our model pre- 
dicts that current experimental parameters can produce 
quadrature squeezing more than 3dB below the zero- 
point motion without requiring an initial ground state or 
perfect detector efficiency. In addition, the measurement 
can be orders of magnitude weaker than that required for 
quantum squeezing via back-action evasion methods. 

The system is described by the Hamiltonian for an os- 
cillator where the spring constant feg = rnuj^ is modu- 
lated with amplitude kr at frequency ujd — ^{ujm + A), 

H=^ + ^[ko + KcosiLJdt + 29)], (1) 

where ujm is the mechanical resonance and m is the effec- 
tive mass. If the continuous measurement of position x is 
fed into a phase-sensitive mixer such as a lock-in ampli- 
fier (see Fig.[T^), as typically done in experiments [l3|, the 
quadrature amplitudes X and Y become the variables of 
interest, where ^Jmujral'hx = X cosiyjmX) + Y "svciiuJmX)- 
The phase Q between the parametric drive and lock-in de- 
tector defines the amplification axis in X-Y space. The 
corresponding quantum operators are X — (a -{- a^)/y2 
and y = — i(a— a^)/\/2. In these variables, Vx = {X"^) = 
Vy = (Y"^) = 1/2 in the ground state due to the un- 
certainty principle. We assume low damping 7 <^ LOm 
corresponding to a high mechanical quality, and small 
perturbation of the spring constant kr <C kg, so that a 
rotating wave approximation can be made in the inter- 
action picture at the frequency ujd, giving 



H = -hAa''a 



2^ 



218-2 

a 



e-^^'o)^) 



(2) 



The second term in this Hamiltonian takes the form of a 
squeezing operator proportional to x = '-^mkr/'^ko. The 
nonlinearity x can be interpreted as the mechanical fre- 
quency shift due to a change in spring constant by kr- 

The conditional evolution of the system, includ- 
ing back-action and thermal noise, is described by 
a stochastic master equation for continuous position 
measurement |22{ . If the measurement rate /i is small 
compared to oj™, the rotating wave approximation can 
be performed on this master equation using the method 
of Ref. [23|. The expectation value of an arbitrary ob- 
servable A is then found to evolve as 
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i2j{N + 1) + fi]{V[a]A) dt (3) 

-^Tfji{n[x]A) dWx + Vw{n[Y]A) dWy , 



where N is the mean thermal phonon number, 77 is 
the quantum efficiency and dWx and dWy are uncor- 
related Wiener processes corresponding to the two po- 
sition quadratures 23| . It is assumed that the measure- 
ment signal also has no thermal fluctuations, which is 
valid for optical readout or a well-cooled signal amplifier. 
The superoperators V and "H arc defined as 

I?[a]i = a'^Aa - -{a^'aA + Aa^'a) 

nialA = aA + Aa^ - (a + a^>i . 

Setting 6* = A = and substituting X and Y for A in Eq. 
([3|) results in damping rates of 7-l-x and 7— x respectively 
for the two quadratures. When there is no nonlinearity 
ix = 0); these reduce to the bare damping rate 7 as 
expected. There exists a threshold for self-oscillation at 
X = 7, where the damping of Y is zero, corresponding 
to infinite gain. At this threshold, the damping of X is 
double its original value, resulting in a gain of 1/2 and 
therefore a noise squeezing limit of 3dB. 

Non-zero detuning A of the pump (and lock-in ref- 
erence) from resonance rotates the phase of the oscil- 
lations, creating elliptical trajectories of the mean val- 
ues and resulting in a shift in the steady-state squeezing 
axis[2J|. Most importantly, as we will show, correlations 
between the quadratures allow the squeezed quadrature 
to be more efficiently localized by measurement. The 
modified below-threshold condition is x^ < A^+7^ where 
above this, the detuning and damping are insufficient to 
keep the mean values decaying to zero, causing instabil- 
ity. With sufficient detuning, this theory is limited only 
by the rotating wave approximation, which requires the 
condition x ^ ^m- 

Near threshold, the squeezing is maximized at an an- 
gle approaching — 7r/4 [2J|. Therefore, for simplicity we 
set the drive phase 6 — n/A so that the squeezing axis 
approximately corresponds to X. The evolution of the 
variances Vx,Vy and covariance C can be derived from 
the master equation using Ito calculus and the relations 
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dC = -d{XY+YX)-(X)d{Y)-{Y)d(X)-d{X)d{Y) 

All terms involving V in Eq. ([3]) , including thermal noise 
and a back-action term /i, add to these position variances. 
The superoperator H, however, results in the conditional 
variance decreasing as the effective measurement rate rjfi 
increases. Setting the differentials above to zero yields 
the steady state solutions as three coupled equations 
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Here, the conditioning parameter z = 8riii'y{N + Nba + 
1/2) characterises how well the motion is detected and 
Nba — /^/27 is the additional phonon number due 
to back-action for continuous measurement. When the 
pump is off, the covariance vanishes and both quadra- 
tures have the standard back-action limited variance con- 
ditioned on the measurement result 
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which in the weak measurement limit (/x ^ 7) gives the 
expected thermal state variance Vt = N + 1/2. The con- 
ditional variance decreases with stronger measurement 
and saturates at the quantum limit of 1/2 in the limit 
/i ^ 7, where the ground state motion can be resolved. 

The parametric pump creates a non-zero covariance, 
or squeezing at an angle of 7r/4 to the X and Y quadra- 
tures, while the detuning rotates the squeezing axis so 
that Vx < Vy- For given parameters, there exists a 
detuning Aopt that optimally uses the measurement to 
minimise Vx- When the measurement is switched off 
{z — 0), this is exactly solvable as Aopt ^ X — 1 with a 
variance limit of Vo/2, agreeing with time domain anal- 
ysis of intracavity parametric squeezing in optics [ITl]. 

A more general analytic result can be obtained for 
the optimal squeezing in the limit that the conditioning 
parameter z 3> 7^, corresponding to high temperature 
and/or strong measurement. Solving for the minimum 
squeezed variance Vx over all detunings results in the 
solution 
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where dashed parameters are normalized by ^77^ + z and 
the optimal detuning Aopt is given by 
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(10) 



It can be seen from Eqs (l51fTU| that A' ^ — x' approaches 
1 in the weak pump limit x ^ 7i giving Vxopt = Vb as 
expected. In the strong pump limit, Aq^j — x' — ^ l/"\/3 
and hence Vx/Vb — > 0. This shows that with quantum 
measurement, improvement on the conditional variance is 
not bounded by 1/2 and perfect squeezing is possible with 
sufficient nonlinearity x'- Fig. [2] compares the analytic 
solutions for Vxopt and Aopt to numerical solutions for 
a fixed measurement strength of ^/7 = 0.1 and various 
initial phonon numbers A'^. The analytic solutions fit well 
at high temperatures and provide a lower bound to the 
squeezing. 
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FIG. 2. Optimal X quadrature squeezing and (inset) opti- 
mal detuning as a function of nonlinearity x ■ Dotted lines: 
numerical solutions for measurement strength /i/7 = 0.1 and 
mean phonon number N = 0,1, 10, 100, 1000. Solid lines: an- 
alytic solutions for no measurement (upper) and for the limit 
.z 3> 7^ (lower). 



The normalized nonlinearity x' is a useful figure of 
merit for classical squeezing at high temperature. In 
the limit where back-action is negligible this can be ex- 
pressed, in terms of mechanical properties, as 



4:y/2NTuipf ko ' 



(11) 



where Q is the mechanical quality factor ujm/j- Hence for 
mechanical squeezing, a high Q and low spring constant 
ko is desirable. If Nrjn/^ is large (i.e. the thermal noise 
is well transduced) the statistics are dominated by the 
measurement conditioning and the effect of the paramet- 
ric drive is reduced. However, lower temperature always 
results in lower absolute variance Vxopt as expected. 

For low initial temperatures and weak measurement, 
z < 7^ so the above solutions are no longer applica- 
ble. Therefore, in order to investigate quantum squeez- 
ing {Vx < 1/2), Eqs (Ulini) were solved numerically for 
initial temperatures near the ground state. The abso- 
lute variance Vx is shown in Fig. [3] for x/7 = 50, com- 
pared to that achievable by BAE (given by Eq. (O where 
Nba = lZI). For relatively weak measurements, the 
DMPA method produces variances well below the zero- 
point motion. This method is also robust to heating, 
with squeezing possible from A^ < 5 at ^/j « 1. At 
this measurement strength, a BAE scheme would require 
A^ < 0.5. Furthermore, the need for an often impractical 
modulation of the measurement strength is eliminated. 

It should be noted that in the limit of very low tem- 
peratures and strong measurement, where back-action 
degrades the DMPA squeezing, BAE is preferable. The 
turning point occurs at a measurement strength of ^/j ~ 
0.4 for A^ = 0. At higher temperatures, stronger mea- 
surement is required to reach the back-action dominated 
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FIG. 3. Optimal quantum squeezing generated by DMPA 
with x/7 = 50 (a,c) and by BAE (b,d). Vxopt is shown as a 
function of fi/y and A'' in (a,b) for J? = 1, and as a function of 
jj,/y and -q in (c,d) for A'^ = 0. Shaded areas denote squeezing 
below the zero-point motion (light) and the standard 3dB 
MPA limit (dark). 



regime. For A^ = and x/l = 50, 6dB of squeezing 
is achievable using DMPA, twice that achievable without 
detuning and measurement. This is robust against detec- 
tion inefficiency, with a reduction to 40% efficiency only 
degrading the squeezing to 5dB. Unconditional squeezing 
can be achieved by using a linear feedback force (e.g. from 
a separate electrode) to stabilize the mean values [2l[. 
This relies on high mechanical Q so that the delay in con- 
verting momentum feedback into physical displacement 
can be neglected. 

We have shown that a detuned parametric amplifier 
with weak measurement can vastly reduce the uncer- 
tainty of one quadrature of motion of a mechanical os- 
cillator. At low initial temperatures, this manifests as 
quantum squeezing. This method therefore opens up the 
possibility of generating mechanical squeezed states in 
many systems, namely those that can be cooled to low 
phonon occupations and parametrically driven indepen- 
dently of the measurement apparatus. 

High frequency mechanical oscillators approaching the 
GHz regime have the well-known advantage of being able 
to be cooled to very low phonon numbers[4]. In addi- 
tion, the small scale involved facilitates a lower intrinsic 
spring constant fco, especially with extreme dimensional 
ratios such those found in oscillators based on carbon 
nanotubes[25|. Many current NEMS resonators are al- 
ready capable of above-threshold parametric amplifica- 
tion and are at high enough frequency to be cooled near 
the ground state. For example, the 9MHz oscillator used 
in Ref. [14| can reach x/7 = 50 with a drive amplitude of 
0.6V. Other suitable examples could involve using elec- 
trostatic forces to drive NEMS resonators coupled to op- 
tical cavities [2, lUl, or even to drive the internal mechan- 



ical modes of optical cavities 26|, |27[ . These setups allow 
resolved sideband cooling of the mechanical mode[28| as 
well as the necessary parametric drive. 
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